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Abstract 

The purpose of this paper is to quantify the interactions of some groups of spoilage organisms that can be usually found in 
refrigerated meat stored in air, such as: Enterobacteriaceae, Pseudomonas, Acinetobacter, Psychrobacter, Shewanella, 
Carnobacterium, Lactobacillus, Leuconostoc, Brochothrix and Kurthia spp. The growth of these organisms was studied in 
the range of temperature 2-1 1°C and pH 5.2-6.4, which is characteristic of refrigerated meat. The main growth parameters 
(maximum specific growth rate and lag time) were modelled by multivariate quadratic polynomials of temperature and pH. 
The interactions of the organisms were analyzed by comparing their growth models obtained in isolation with those obtained 
in mixture. The difference between the models was quantified by statistical F -values which were used to measure how much 
the growth of an organism or group of organisms was affected by others and which of them dominated their joint growth. 
© 1998 Elsevier Science B.V. 

Keywords: Predictive microbiology; Spoilage; Interactions; Pseudomonas; F-test 



1. Introduction 

Predictive food microbiology has been focusing on 
modelling the microbial responses to the food en- 
vironment, in the interest of food safety and of 
avoiding spoilage. Predictive models have been 
regularly published to describe the growth of a strain 
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or a mixture of strains as a function of the en- 
vironmental conditions characteristic of food. 

Because of the variety and number of spoilage 
organisms, spoilage models are less straightforward 
to develop than pathogen models and their applica- 
tion is much more limited. Dalgaard (1995) sug- 
gested that predictive models for spoilage should be 
developed only after knowing the micro-organisms 
responsible for the reactions important in the process 
of spoilage and the range of environmental con- 
ditions under which these organisms cause spoilage. 
McMeekin and Ross (1996) emphasised that spoil- 
age models, based on the responses of the dominant 
organism, are valid only in a specific range of 
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conditions. Out of this range, different bacteria or 
metabolites may be responsible for the spoilage, 
making the model no longer valid. 

The process of spoilage in aerobicaJly stored 
chilled meat is well characterised and this food 
system was chosen for our work. Pseudomonas, 
Acinetobacter and Psychrobacter spp. are the domi- 
nant bacteria causing spoilage, when their concen- 
tration reaches - 10 7 /cm 2 , or /g, as a consequence 
of their fast growth and their metabolism. Brochoi- 
hrix thermosphacta, cold- tolerant Enterobac- 
teriaceaea and lactic acid bacteria also occur in 
aerobically-stored meat but generally show slower 
growth and /or cause less offensive or later spoilage 
signs (Gill and Newton, 1977; Dainty and Mackey, 
1992; Lambropoulou et ah, 1996). 

The growth conditions that can be found in chilled 
meat are: high water activity (-0.99); pH values of 
-5.8 though they can go from -5.4 to -7.0, and 
refrigeration temperatures (Egan and Roberts, 1987). 
The process of chilled meat spoilage can be iden- 
tified as the growth of the aerobic gram-negative 
bacteria or simply of Pseudomonas spp. 

Most foods are complex systems with heteroge- 
neous populations. A mechanistic approach describ- 
ing the interactions of different strains in food would 
require a complex mathematical model with many 
variables. Here, we do not make an attempt to 
describe the dynamic coexistence of organisms in 
meat but we take a simpler approach. We analyse the 
differences between the isolated and joint growth of 
four groups of spoilage organisms, that can be 
usually found in refrigerated meat stored aerobically. 
The results enables us to quantify how much the 
growth of a group is affected by the other groups and 
to what extent one group dominates the spoilage 
flora. We define a subpopulation of a mixture of 
organisms dominant if its specific growth rate is 
higher than that of the other subpopulations. A 
dominant subpopulation does not always outnumber 
the others because the actual number of cells de- 
pends also on the lag phase and the inoculum size. 
However, these latter parameters depend on the 
initial conditions (the lag depends on the initial 
physiological state) and, in fact, the maximum 
specific growth rate is the only intrinsic parameter of 
bacteria] growth. If one of the subpopulations has 
higher maximum specific growth rate than the others, 
then it outnumbers the others provided approximate- 



ly equal inoculum size and initial physiological state 
at the start. 

After comparing the growth models, a generic 
spoilage model is created for the specific growth 
rate, by using the data of the groups found as 
dominant. In accord with the above note on the 
concept of dominance, the generic model can help to 
predict the time to spoilage assuming that the 
dominant organisms are present in sufficient number 
at all. 



2. Materials and methods 

2.1. Bacterial strains 

Thirty-two strains of Pseudomonas, Acinetobacter, 
Psychrobacter, Shewanella, Enterobacteriaceae, Car- 
nobacterium, Lactobacillus, Leuconostoc, B. ther- 
mosphacta and Kurthia spp., were used in this work. 
They were divided into four groups as follows. 

Group P: Pseudomonas putida NCFB 754 (spoiled 
milk), Pseudomonas fragi NCFB 2902 (beef), Pseu- 
domonas lundensis NCFB 2908 (minced beef), 
Shewanella putrefaciens NCFB 756 (tinned butter), 
Acinetobacter sp. NCIMB 11168 (lamb carcass 
meat), Acinetobacter sp. NCIMB 11 169 (lamb car- 
cass meat), Psychrobacter inmobilis NCIMB 11372 
(pork sausage), Psychrobacter immobilis NCIMB 
11650 (meat). 

Group E: Enterobacter agglomerants NCFB 2071 
(minced meat), Enterobacter agglomerants NCFB 
2073 (pasteurised milk), Klebsiella oxytoca NCFB 
2678 (Cheddar cheese), Klebsiella pneumoniae 
NCFB 1010 (unknown), Escherichia coli NCFB 555 
(raw milk), Proteus morganii NCTC 235 (unknown). 

Group L: Lactobacillus sp. NCFB 2812 (vacuum 
packed pork), Lactobacillus sp. NCFB 2813 (vacuum 
packed beef), Lactobacillus sp. NCFB 2814 (vacuum 
packed bacon), Lactobacillus sp. NCFB 2815 (vac- 
uum packed lamb), Leuconostoc carnosum NCFB 
2775 (vacuum packed meat), Leuconostoc gelidum 
NCFB 2776 (vacuum packed meat), Leuconostoc 
gelidum NCFB 2800 (vacuum packed beef), Car- 
nobacterium divergens NCFB 2856 (vacuum packed 
beef), Carnobacterium divergens NCFB 2857 (vac- 
uum packed lamb), Carnobacterium piscicola NCFB 
2853 (vacuum packed beef), Carnobacterium pis- 
cicola NCFB 2854 (vacuum packed lamb). 



C. Pin, J. Barony i I International Journal of Food Microbiology 41 (1998) 59-72 



61 



Group B: Brochothrix thermosphacta NCFB 1676 
(fresh pork sausage), Brochothrix thermosphacta 
NCFB 2830 (lamb), Brochothrix thermosphacta 
NCFB 2849 (vacuum packed pork), Kurthia gibsonii 
NCIMB 10499 (pork sausage), Kurthia zopfii 
NCIMB 10494 (hamburgers), Kurthia zopfii NCIMB 

10496 (frozen mince pork), Kurthia zopfii NCBvlB 
10498 (pork sausage), Kurthia gibsonii NCIMB 
10495 (fat trimmings), Kurthia gibsonii NCIMB 

10497 (rendered lard). 

2.2. Media 

Tryptone Soya Broih (TSB, Oxoid/Unipath 
CM 129) was made in 3000-ml volumes, adjusted to 
the target pH values using 5 M HC1, dispensed in 
250-ml volumes and autoclaved at 121°C for 15 min. 

The bacterial counts were obtained on: Tryptone 
Soya Agar for non-selective counts (TSA, Oxoid CM 
131); Cetrimide Fusidin Cephaloridine agar for 
Pseudomonas spp. (CFC, Oxoid CM559, SR 103); 
MacConkey agar for Enterobacteriaceae (Oxoid/Un- 
ipath CM115); MRS agar for lactic acid bacteria 
(MRSA, Oxoid CM 361); Streptomycin Thallous 
Acetate Actidione agar for B. thermosphacta (STAA, 
Oxoid CM 881, SR151). 

2.3. Inoculum preparation 

The cultures of each strain were activated by 
transferring loop inocula grown in 10 ml of TSB at 
25°C at three successive 24-h intervals immediately 



prior to the experiments. Equal volumes of the 
cultures, prediluted approximately to 10 cfu/ml, 
were combined to give cocktails of the groups P, E, 
L, and B. Equal volumes of this cocktail were used 
to get a total mixture (T) containing all the strains 
together. 

2.4. Experimental procedure 

TSB was maintained overnight at the intended 
storage temperature. After measuring the pH, the 
broth was inoculated with the appropriate cocktail 
giving a final concentration of ~ 10 3 cfu/ml. Imme- 
diately after inoculation, the broth was aseptically 
dispensed in 10-ml volumes, using a peristaltic pump 
(Accuramatic Mk 5; Jencons Ltd., Leighton Buzzard, 
UK), into sterile 1-oz universal bottles with plastic 
screw caps. The initial inoculum level was deter- 
mined from one of the bottles and the remaining 
bottles were incubated at the desired temperature 
(±6.5°C) in static conditions. During the storage, the 
samples were plated onto either selective or non- 
selective agars, according to Table 1. Plates were 
incubated for 24-48 h at 25°C, using anaerobic jars 
(Oxoid/Unipath, HP11) with Gas Generating Kits 
(Oxoid/Unipath BR038B) when required. 

Growth responses were studied under 12 con- 
ditions, combining four temperatures (2, 5, 8 and 
11°C) and three pH values (5.2, 5.8 and 6.4), as 
shown in Table 1. 

The ability of a selective agar to recover a certain 
strain was checked by growing the strain at 25°C, for 



Org^Lon of the 'type A', 'type B' model, and the total model. Each mode! is based on growth curves generated at the 12 combinations 
of four temperatures (2. 5, 8, \ 1°C) and three pH values (5.2, 5.8 and 6.4) 



Model 



Group of organisms 



Growth media 
in broth 



Plating 
media 



Plates 
incubation 



Type A 


couc 

M(P) 
M(E) 
M(L) 
M(B) 


P 
E 
L 
B 


TSB 
TSB 
TSB 
TSB 


TSA 
TSA 
MRSA 
TSA 


Aerobic 
Aerobic 
Anaerobic 
Aerobic 


Type B 


M(T,CFC) 
M(T.MacC) 
M(T,MRSA) 
M(T,STAA) 


TotaJ (all the four groups inoculated together) 


TSB 
TSB 
TSB 
TSB 


CFC 

MacConkey" 

MRSA 

STAA 


Aerobic 
Aerobic 
Anaerobic 
Aerobic 


Total model 


M(T.TSA) 


Total 


TSB 


TSA 


Aerobic 



'Only red colonies were counted to select group E. 
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24 h, then plating it on the selective agar (Table 2, 
first five columns). 

Each strain was incubated at the four temperatures 
for 30 days to establish whether they were able to 
grow at those temperatures at all (Table 2, last 
column). 

2.5. Growth models 

The generated bacterial growth curves were fitted 
by the model of Baranyi and Roberts (1994) (see 
Fig. la and b). The main growth parameters, esti- 
mated at each combination of temperature and pH, 
were: 

• maximum specific growth rate (p.) 

• lag phase (X). 

When a mixture of organisms grow together, the 
specific growth rate is not constant, even if all the 
subpopulations are in the exponential phase. The 
instantaneous specific growth rate of the i-th sub- 
population, at the time /, is 



where *.(f) is the respective bacterial concentration. 
Let the overall concentration be denoted by 

*(0 = *,(0+* 2 «)+ ... 

Then, as it can be derived mathematically, the 
instantaneous specific growth rate of the whole 
population, at the time t is: 

i.e. the weighted average of the individual specific 
rates, where the weights are the actual proportions of 
the individual sub-concentrations in the whole popu- 
lation. Therefore, assuming that the fastest growing 
subpopulation does not have longer lag and smaller 
starting number than the others, the dominance in 
rate means numerical dominance very soon and the 
specific rate of the whole population becomes practi- 
cally indistinguishable from the fastest specific 
growth rale. This justifies the use of the model of 
Baranyi and Roberts (1994), to fit growth curves of 
mixed cultures, since that model is based on the 



assumption that the specific growth rate is practically 
constant for a phase. 

To correct the heterogeneity of the variance, the 
natural logarithm (In) of u- and X was modelled as a 
function of temperature and pH. This model was a 
standard multivariate second order polynomial as 
described by McClure et al. (1993). About the use of 
the variance-damping Iog-transformauon, see Alber 
and Schaffner (1992). 

Each of the nine growth models created was 
denoted by a code referring to the name of the 
modelled group and possibly to the used selective 
media (Table 1). The model codes do not include 
any reference to the modelled growth parameter and, 
unless we explicitly state in another way, the term 
'model' refers to In p.. 

Two types of predictive models characterising the 
four groups and a total model were generated for 
comparisons. 

Type A: a group inoculated in isolation. The four 
models of 'type A' were denoted by M(P), M(E), 
M(L) and M(B), respectively. 
Type B: all the above groups inoculated together 
in a total mixture (T) and plated on selective agar. 
The four models of 'type B* were denoted by 
M(T,CFC), M(T,MacC), M(T,MRSA) and M(T- 
,STAA), respectively. 

Total model: all the above groups inoculated 
together in a total mixture (T) and plated on 
non-selective, ordinary TSA. The model obtained 
this way is denoted by M(T.TSA). 

After comparing the growth models, also a generic 
spoilage model, M(G), was created for the specific 
growth rate, by using the data of the groups found as 
dominant. 

2.6. Measuring microbial interactions 

(a) Differences between the growth rates of the 
groups in isolation were studied by comparing the 
models M(P), M(E), M(L), and M(B) with each 
other. 

(b) To study how the growth of a particular group 
was affected by the total mixture, its polynomial 
model for In p. ('type A') was compared with the 
model obtained on its respective selective medium 
from the total mixture ('type B'): 
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M(P) vs. M(T,CFC) 

M(E) vs. M(T,MacC) 

M(L) vs. M(T,MRSA) 

M(B) vs. M(T,STAA) 

(c) To identify the dominant group in the total 
mixture, the models for both In \l and X obtained on 
selective media ( 4 type B') were compared with those 
obtained on the non-selective TSA (Total model): 

M(T,CFC) 
M(T, TSA) vs. M(T,MacC) 
M(T,MRSA) 
M(T,STAA) 

2.7. F test 

An F test was used to decide whether the differ- 
ence between two models was significant. The 
procedure of the test can be summarised as follows. 

Suppose that a series of n observations, y, . . . y„, 
(in our case, growth rates or lag times under various 
environmental factors) are fitted by the model M l 
with p { parameters. Similarly, let a series of m 
observations, y n + x .-.y„+ m , be fitted by the model 
M 2 with p 2 parameters, and let M 2 have similar 
structure (such as polynomial) as M,. Furthermore, 
let the merged n+/n observations, y, . . . y„ +m , be 
fitted by a unified model, M M , with p u parameters, 
and with the same model structure as M, and M 2 . In 
other words, wc have 

1. a separated model (M s ) of p s =P, +P 2 parameters 
containing two segments, M, and M 2 , laid over 
the observations >>, . . . y n and y n+I ... y n+m , re- 
spectively; 

2. a unified model (M u ) of p u parameters laid over 
the unified observation points, y, . . . y H + m - 

Under general conditions imposed on the variance 
of the data, the question, whether the difference 
between the models M, and M 2 is significant, can be 
tested by the F value (Box and Draper, 1987) 

X(y.-y.) 2 . 
2<y,-'fl) J 

/"I 

n + m- p u 



where y, is the fitted value for the /-th observation 
when the separated model, M,, is used for fitting and 
y. is the fitted value when the unified model, M w , is 
applied. 

When the probability belonging to an F value was 
less than 0.05 (P<0.05), the difference between the 
two models was considered significant. 



3. Results 

3.1. Recovery of the organisms on selective media 

The recovery of the different organisms on the 
used selective media is shown in Table 2. 

CFC was used to recover the group P from the 
total mixture. The Acinetobacter and Psychrobacter 
strains of the group P did not grow on CFC but did 
grow on MacConkey agar giving red colonies. Three 
strains from the group E were able to grow on CFC 
agar. Among the bacteria of the group P, the 
Pseudomonas and Shewanella strains grew on CFC, 
while the Acinetobacter and Psychrobacter strains 
did not. The group P grown alone and the total 
mixture plated on CFC showed similar results. This 
indicates that the Pseudomonas and Shewanella 
strains were primarily responsible for the growth 
within this group when inoculated separately. This 
was supported by the fact that one Acinetobacter 
strain grew neither at 2°C nor at 5°C and one 
Psychrobacter strain did not grow at 2°C (Table 2). 
Hence the failure to recover the Acinetobacter and 
Psychrobacter strains on CFC did not affect the 
comparison between M(P) and M(T,CFC). 

Moreover, the Acinetobacter and Psychrobacter 
strains gave red colonies on MacConkey agar. As it 
will be shown, the growth of those bacteria giving 
red colonies on MacConkey agar was thoroughly 
inhibited in the total mixture. Because of this strong 
inhibition, the mentioned four strains of the group E 
did not have any effect on the counts obtained on 
CFC agar. 

For these reasons, Pseudomonas and Shewanella 
strains were considered to be the organisms respon- 
sible for the growth of the group P when inoculated 
alone. Besides, CFC agar could be used to measure 
the growth of the group P in the total mixture. 

MacConkey agar was used to recover the group E 
from the total mixture, counting red colonies only. 
Pseudomonas and Shewanella strains (from the 
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Tabic 2 . . 

Recovery of the different strains of the four groups on selective media 



Microorganisms 



"Logcounts of culture, grown at 25°C for 24 h, on: Lowest growth temperature studied CQ 



Group P: 

Pseudomonas purida NCFB 754 
Pseudomonas fragi NCFB 2902 
Pseudomonas lundensis NCFB 2908 
Shewanella putrefaciens NCFB 756 
Acinetobacter sp. NC1MB 1 1 168 
Acinetobacter sp. NC1MB 11 ! 69 
Psychrobacter inmobilis NCIMB 1 1372 
Psychrobacter inmobilis NCIMB 1 1650 

Group E: 

Enterobacter agglomerants NCFB 207 1 
Enterobacter agglomerants NCFB 2073 
Klebsiella oxytoca NCFB 2678 
Klebsiella pneumoniae NCFB 1010 
Escherichia coli NCFB 555 
Proteus morganii NCTC 235 

Group L: 

Camobacterium divergens NCFB 2856 
Carnobacterium divergens NCFB 2857 
Camobacterium piscicola NCFB 2853 
Camobacterium piscicola NCFB 2854 
Lactobacillus sp. NCFB 2812 
Lactobacillus sp. NCFB 2813 
Lactobacillus sp. NCFB 2814 
Lactobacillus sp. NCFB 2815 
Leuconostoc camosum NCFB 2775 
Leuconostoc gelidum NCFB 2776 
Leuconostoc gelidum NCFB 2800 

Group B: 

Kurthia gibsonii NCIMB 10499 
Kurthia zopfii NCIMB 10494 
Kurthia zopfii NCIMB 10496 
Kurthia zopfii NCIMB 10498 
Kurthia gibsonii NCIMB 10495 
Kurthia gibsonii NCIMB 10497 
Brochothrix thermosphacta NCFB 1676 
Brochothrix thermosphacta NCFB 2830 
Brochothrix thermosphacta NCFB 2849 



TSA 


CFC 


McConkey 


MRSA 


STAA 




9.6 


9.5 


9.5 (R + Y) 


ND 


ND 




10 


10 


9.7(Y) 


ND 


ND 


2 


10 


9.9 


9.9(Y) 


INU 


ND 
ySXJ 


2 


9.8 


9.9 


9.9(Y) 


Kin 

INU 




2 


9.6 


9.0 


ND 


Kir* 




2 


9.8 


ND 


9.6(R) 


INU 


INU 


g 


9.8 


ND 


9.8(R) 


NU 


Mn 

INLJ 


2 


9.7 


ND 


9.7(R) 


ND 


Kin 

INU 


9 
i. 


9.6 


ND 


9.6XR) 


ND 


Kin 
INU 




10 


9.0 


9.9(R) 


5.25 


ND 




9.7 


ND 


9.5(R) 


6.28 


ND 


t 


9.4 


ND 


9.5(R) 


ND 


ND 


z 


in 

1U 


9.9 


10(R) 


ND 


ND 


2 


10 


10 


10(R) 


ND 


ND 


2 


10 


ND 


9.9(R) 


ND 


ND 


Q 

o 


9.4 


9.7 


9.4{Y) 


ND 


ND 


2 


9.1 


ND 


ND 


9.3 


ND 




8.8 


ND 


ND 


9.2 


ND 


I . 


9.3 


ND 


ND 


9.4 


ND 


*> 
2 


9.0 


ND 


ND 


9.2 


ND 


2 


9.3 


ND 


ND 


9.5 


ND 


2 


9.2 


ND 


ND 


9.3 


ND 


L 


9.6 


ND 


ND 


9.6 


Nu 


L 


9.1 


ND 


ND 


9.3 


ND 


2 


9.3 


ND 


ND 


9.4 


ND 


5 


9.2 


ND 


ND 


9.5 


ND 


2 


9.7 


ND 


ND 


9.4 


ND 


2 


9.3 


ND 


ND 


9.4 


ND 


2 


9.2 


ND 


ND 


ND 


92 




8.4 


ND 


ND 


ND 


ND 


5 


8.3 


ND 


ND 


ND 


ND 


2 


9.0 


ND 


ND 


ND 


ND 


2 


8.7 


ND 


ND 


ND 


ND 


2 


8.7 


ND 


ND 


ND 


ND 


2 


9.2 


ND 


ND 


ND 


ND 


2 


8.8 


ND 


ND 


ND 


8.7 


2 


8.3 


ND 


ND 


ND 


8.0 


2 


8.7 


ND 


ND 


ND 


9.0 


2 



ND, growth not detected; (Y), yellow colonies; (R), red colonies. 



group P) gave yellow colonies on MacConkey agar. 
This made it difficult and inaccurate to count the red 
colonies. In spite of all these problems, the results on 
MacConkey agar let us measure the growth of the 
group E in the total mixture. 

MRSA worked quite well as a selective medium 
for the group L. 



On STAA, Kurthia spp. could not grow unlike B. 
thermosphacta. When the group B was grown in 
isolation, Kurthias spp. did not grow. This was 
concluded from the result that we observed the same 
growth curves when incubating the plates of the 
group B under aerobic (both Brochothrix and Kurth- 
spp. were able to grow) and anaerobic (only 



ta 
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Brochothrix spp. could grow) conditions. Hence, the 
growth of the group B was attributed exclusively to 
B. thermosphacta, and STAA was used to measure 
the growth of the group B in the total mixture. 

3.2. Isolated growth of the groups (data for 
models of 'type A ') 

The maximum specific growth rates (pO and the 
lag times (X) of the groups P, E, L and B, when 
grown in isolation, were derived by the curve fitting 
procedure of Baranyi and Roberts (1994). The 
logcount curves of the four groups, grown at 2°C and 
pH 5.8, are shown in Fig. la, for an example. The 
fitted growth parameters, for every studied combina- 
tion of temperature and pH, are shown in Table 3. 

Fig. la suggests that the group L is slower than 
the other three. The group P has the fastest growth 
but its rate is just slightly higher than that of the 
group B which, on the other hand, has the longest 
lag. The question to be answered is whether these 
observations are true for other temperature and pH 
values. 

Fig. 2a represents, in 3-D, the maximum specific 
growth rates of Table 3 as a function of temperature 
and pH. It shows that the observation made about the 



group L in Fig. la, namely that it grows more slowly 
than the others, is true only at low temperatures. In 
the main, however, Table 3 and Fig. 2a confirm what 
Fig. la demonstrated about the difference between 
the growth rates of the four groups. 

When comparing the 'type A' models with each 
other (specific growth rate models obtained in isola- 
tion), the F-test results showed no significant differ- 
ences between group P and E (P=0.33). The specific 
growth rate of group B was also similar. The 
difference between groups B and E was not signifi- 
cant (P = 0.37). Between groups P and B, it was 
significant at 5% but not at 1% (/> = 0.021). The 
model of group L, however, was significantly differ- 
ent (P<0.05) from each of the other three models. 

3.3. Growth of the groups in the total mixture 
(data for models of 'type B ') 

Fig. lb shows the fitted growth curves produced 
by the total mixture when plated on four different 
agars each selective for the four groups. This figure 
also shows a fifth fitted curve which was obtained on 
the non-selective TSA, from the total mixture. 

Fig. lb suggests that the growth of the group P 
(selected, from the total mixture, by CFC agar) is not 





Fig. 1. (a) Fitted logcounts obtained when each of the four groups was grown in isolation, at 2«C pH 5.8 Curves hke these give ^ basis fo 
the mode s of -type A' group P; A. group E; group L; group B. (b) Fitted logcounts obtained when the total nature of the four 
groups was grown at 2°C. pH 5.8 and plated on the four selective agars (thin lines) and on TSA (thick line). Curves like these give the basis 
fa L modds of Type B' (thin lines) and for the total model (thick line). The name of the ^oup for which an agar « ^ ^cuve js in 
parenthesis in the figure. The legend is parallel to (a), according to the group-agar relauon: O, CFC; A, MacConkey; O, MRSA, □, STAA, 
X T total mixture plated on TSA. 
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Table 3 

Maximum specific growth rales (l/h) and 
data for the models of 'type A 



C. Pin, 7- Baranyi I International Journal of Food Microbiology 41 (1998) 59-72 

limes (h) of ihe groups inoculated in isolation. THese growth parameters provided ihe input 



Temp 

rc) 



pH Group P 



Group E 



Group L 



Rcl. 



Ret 

or. 



ReL 



RcL 
cn. 



RcL 
en. 



2 6.4 

2 5.8 
2 

5 M 

5 5.8 

5 5.2 

8 6.4 
8 5.8 
8 5.2 



6.4 
5.8 
51 



0.090 
0.084 
0.076 

0.12 
0.11 
0.11 

0.20 
0.19 
0.15 

0.28 
0.27 
0.22 



7.9* 
5.7% 
6.5% 

110% 
6J% 
8.3% 

7.0% 
7.0% 
3.3% 

7.9% 
5.5% 
8.3% 



30 


22% 


0.071 


4.2% 


37 


110% 


29 


18% 


0.072 


3.0% 


36 


8.0% 


48 


13% 


0.063 


3.7% 


44 


8.4% 


32 


21% 


0.081 


4.3% 


24 


19.0% 


29 


15% 


0.077 


3.9% 


20 


210% 


46 


11% 


0.077 


4.2% 


41 


8.9% 


16 


13.0% 


0.20 


9.6% 


24 


10.0% 


16 


14.0% 


0.17 


12.0% 


24 


110% 


18 


7.9% 


0.17 


7.6% 


28 


7.7% 


7.0 


26% 


0.26 


9.7% 


12 


12.0% 


43 


36% 


0.26 


6.5% 


11 


10.0% 


6.8 


41% 


0.22 


5.6% 


11 


14.0% 



0.044 
0.042 
0.034 

0.066 
0.059 
0.049 

0.14 
0.14 
0.12 

019 
0.29 
0.27 



3.6% 
2.8% 
17% 

3.9% 
17% 
3.7% 

3.3% 
14% 
11% 

6.0% 
41% 
3.1% 



31 
30 
26 

21 

20 
23 

10 
10 
II 



Group I 



err. 



25% 
22% 
33% 

28% 
21% 
33% 

21% 
14% 
13% 



12 95% 
19 50% 
18 43% 



Rel 
err. 



ReL 

err. 



0.11 

0.084 

0.068 

0.11 

0.086 

0.079 

0.18 
0.16 
0.14 

016 
0.24 
0.31 



6.8% 
7.1% 
3.6% 

3.5% 
4.9% 
16% 

19% 
3.0% 
4.9% 

7.3% 
3.6% 
13.0% 



67 


6.3% 


98 


6.2% 


89 


41% 


44 


6.5% 


64 


71% 


110 


13% 


18 


8.5% 


28 


5.7% 


43 


61% 


8.7 


29% 


10 


13% 


27 


10% 



jl, maximum specific growth rate (l/h). 
X, Ug rime (h). 

ReL err., standard error of the fitted 



expressed as percentage of the parameter. 



different from that of the total mixture plated on 
TSA. The other groups grew slower than the group P, 
and the group E did not grow in the experimental 
time at all. 

The growth parameter estimates, as produced by 
the curve fitting described above, can be seen in 



Table 4, for every studied combination of tempera- 
ture and pH. That table, with Fig. 2b, demonstrates 
that the observations made on the specific growth 
rates in Fig. lb are typical for the whole region of 
environmental factors analyzed here. 
While the group P was the dominant in the mixed 





A, MacConkey; O, MRSA; STAA; X, lotal mixture plated on TSA. 
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Table 4 

Maximum specific growth rates (1 /h) and lag times (h) of the total mixture plated on selected agars (first four sections) and on TSA (last 

section). Tfcese growth parameters provided the input data for the models of type B and the total model, respectively 

S1AA(B) TSA(non-selecuve) 



Temp 


pH 


CFC(P) 








MacC(E) 






MRSA(L) 






rc> 




V> 


RcL 


X 


Rel. 




RcL 


X 


Rel 


V- 


Rel. 


x 


RcL 










err. 




err. 




en. 




en. 




en. 




6.4 


0.065 


7.1% 


32 


31% 










0.062 


4.3% 


59 


10.0% 




5.8 


0.087 


7.8% 


57 


10% 










0.060 


3.6% 


66 


7.5% 




5.2 


0.075 


15.0% 


64 


20% 










0.056 


16% 


73 


5,1% 




6.4 


0.10 


6.0% 


24 


17% 










0.071 


4.5% 


32 


16.0% 




5.8 


0.11 


8.4% 


35 


13% 










0.076 


5.6% 


45 


110% 




5.2 


0.080 


9.1% 


29 


30% 










0.072 


3.8% 


43 


10.0% 




6.4 


0.17 


5.1% 


13 


20% 


0.10 


11.0% 


12 


26% 


0.13 


18% 


5.4 


59.0% 




5.8 


0.19 


6.8% 


20 


11% 


0.076 


4.8% 


12 


28% 


0.12 


2.4% 


7.9 


28.0% 


8 


5.2 


0.17 


5.7% 


22 


9.0% 


0.063 


7.9% 


85 


63% 


0.12 


10% 


8.9 


20.0% 


11 


6.4 


0.25 


7.8% 


8.5 


27% 


0.21 


4.3% 


9.3 


15% 


0.20 


2.2% 


4.7 


27% 


11 


5.8 


0.27 


8.6% 


9.4 


21% 


0.20 


5.3% 


7.0 


30% 


0.21 


2.2% 


4.2 


30% 


11 


52 


0.25 


6.3% 


10 


15% 


0.18 


4.4% 


8.1 


20% 


0.20 


2.0% 


4.1 


29% 



RcL X Rel. ji Rel. X Rel 
en. err. err. err. 



0.072 4.7% 50 12% 0.057 5.5% 23 44% 

0.058 9.4% 62 16% 0.086 10.0% 66 10% 

0.040 15.0% 82 16% 0.074 11.0% 74 13% 

0.085 3.9% 19 28% 0.089 4.2% 21 15% 

0.069 7.3% 23 33% 0.097 6.5% 34 11% 

0.055 19.0% 43 30% 0.075 8.4% 31 28% 

0.18 4.0% 12 16% 0.16 5.1% 12 22% 

0.16 5.1% 12 22% 0.17 5.5% 17 14% 

0.10 5.6% 8J 46% 0.15 6.5% 18 15% 

0.26 4.1% 5.7 28% 0.27 10.0% 9.8 21% 

015 10.0% 8.6 34% 0.25 8.2% 85 25% 

011 12.0% 8.9 34% 0.24 9.4% 9.3 25% 



population for all of the conditions, the group E did 
not grow at 2°C and 5°C. 

3.4, Comparing models of 'type A \ 'type B* and 
the total model 

To measure how other organisms affected the 
growth of a particular group, its maximum specific 
growth rate was compared with those obtained on 
selective agar from the total mixture. This com- 
parison was carried out on the specific growth rate 
values associated to the filled and hollow symbols in 
Fig. 2a and b, respectively. 

To measure which group was the dominant in the 
total mixture, the maximum specific growth rates 
obtained on the respective selective media are to be 
compared with those obtained on TSA from the total 
mixture (see the difference between die hollow 
symbols and the 'cross* symbol in Fig. 2b). Lag 
phases were analyzed in the same way. 

The comparison is even more efficient when 
carried out on fitted response surfaces. Those were 
generated from the data represented by: 

— the hollow symbols of Fig. 2a 

(gaining 'type A* models: M(P), M(E), M(L) 
and M(B)); 

— the filled symbols of Fig. 2b 



(gaining 'type B' models: M(T,CFC), 
M(T,MacC), M(T,MRSA) and M(T,STAA)); 
— the 'cross' symbols of Fig. 2b 

(gaining the total model: M(T,TSA)). 

The fitted coefficients of the above quadratic 
models for In \l and In X, as well as the respective 
standard errors of fitting are tabulated in Table 5. 

To establish whether the difference between two 
response surfaces was significant, the F test was 
applied. By that, for each group, we could obtain 
quantified answers lo the following questions. 

1. Did the group grow in a significantly different 
way in the total mixture than in isolation? (Com- 
paring its 'type A' model with the respective, 
'type B' model.) 

2. Can its growth be identified with the total 
growth? (Comparing its 4 type B' model with the 
total model.) 

The F- value can be considered as a certain 
'distance* between two models, therefore the higher 
the F-value is, the more reason we have to believe 
that two compared models are different. High F~ 
values indicate negative answers to the above two 
questions, which means, in our case: (i) inhibitory 
effect of the other organisms on the group in 
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Efficients of the quadratic surfaces fined co the natural logarithm of the specific growth rates (In rt ™* ** ^ (In X) 



Model 
code 



Growth Coefficients 



parameter Const 



Temp 



pH 



Temp.pH 



Temp" 



PH 2 



S.E of 
fitting 



Type A M(P) 
M(E) 
M(L) 
M(B) 

Type B M(T,CFC) 
M(T,MacC) 
M(T,MRSA) 
M(T.STAA) 

Total model M(T,TSA) 

M(G)* 



Ln p. 
Ln X. 
Ln p. 
Ln X 
Ln u. 
Ln X 
Ln u. 
Ln X 

Ln m- 
Ln X 
Ln u. 
Ln X 
Ln p. 
Ln X 
Ln p. 
Ln X 

Ln p. 
Ln X 

Ln p. 



— o.yi i 




1.369 


0.01182 


31.39 


-0.1110 


-9.231 


0.05013 


-4.777 


0.04475 


0.5848 


0.008529 




-0.1508 


-5.830 


0.03131 


-9.409 


0.2006 


1.780 


-0.01923 


-1.469 


0.3519 


1.428 


-0.03801 


0.1956 


0.2776 


-1.476 


-0.05605 


10.84 


0.3834 


-1.974 


-0.07063 


-16.90 


0.0633 


4.920 


0.005542 


-8.734 


-0.3371 


5.047 


0.03325 


-5.833 


0.3708 


0.02260 




2.097 


-0.07262 


0.1365 




-6.184 


0.06183 


1.056 


-0.004472 


0,1225 


-0.4628 


2.021 


0.02395 


-14.08 


0.2029 


3.252 


-0.02444 


4.847 


-0.6598 


0.5918 


0.03870 


-14.14 


-0.1173 


4.112 


0.02870 


-7.275 


-0.7663 


5.060 


0.09737 


- 12.65 


-0.004356 


3.467 


0.01535 



0.001505 
-0.02932 

0.004685 
-0.01194 

0.01035 
-0.03041 

0.01390 
-0.01438 

0.003724 
-0.002973 



0.008024 
-0.0004022 
0.008459 
0.01485 

0.007475 
0.001154 

0.004234 



-0.1109 
0.7517 

-0.04605 
0.4709 

-0.1295 

-0:1066 
0.1750 
0.1533 

-0.4270 
-0.4827 



-0.08525 
-0.2034 
-0.2324 
-0.09711 

-0.3694 
-0.5249 



0.066 
0.122 
0.209 
0.211 
0.099 
0.074 
0.109 
0.207 

0.133 
0.121 
0.149 
0.140 
0.062 
0.396 
0.134 
0.282 

0.136 
0.179 



-0.3024 0.126 



•Generic spoilage model for In p based on the growth rates of the dominant organisms. 



question; (ii) that the considered group is not among 
the dominant ones in the total growth. 

3.4.1. Results comparing 'type A' and 'type B' 
models 

Fig. 3 and Fig. 4 compare the models M(P) and 
M(E) with M(T,CFC) and M(T,MacC), respectively. 
Note that, in the case of M(T,MacC), the response 
model is a linear surface because growth was 
delected at the highest two temperatures only. From 
Fig. 3 and Fig. 4, it can be concluded that the growth 
of the group P was not affected by the other groups, 
but the group E grew significantly slower in the total 
mixture than in isolation. 

The distances between the models were expressed 
in terms of the F test results, as can be seen in Fig. 
5. 

Fig. 5 demonstrates that the models of the groups 
E t L and B, when grown in isolation (denoted by 
M(E), M(L) and M(B), respectively), showed sig- 
nificantly faster growth than those obtained on 
MacConkey, MRS A and STAA from the total mix- 
ture (denoted by M(T,MacC), M(T,MRSA), and 



M(T,STAA), respectively). This confirms that there 
is no significant difference between the growth 
models of the group P, whether grown in isolation or 
in the total mixture, and that the biggest discrepancy 
was found between the two models of the group E, 




Fig 3. Comparing the growth of group P in isolation with that in 
the total mixture. (Models M(P) and M(T.CFC)). Legend as in 
Fig. la and Fig. lb. Thick surface: M(P). Thin surface: 
M(T.CFC). 
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Fig. 4. Comparing the growth of group E in isolation with that in 
the total mixture. (Models M(E) and M(T.MacC)). Legend as in 
Fig. la and Fig. lb. Thick surface: M(E). Thin surface: 
M(T\MacC). 

indicating that the other groups severely affected the 
growth of group E. 



3.4.2. Results comparing 'type B ' models and the 
total model 

The dominant group responsible for the growth in 
the total mixture was identified by determining not 
only the group or groups with the fastest growth rate 
but also with the shortest lag phase. Therefore, the F 

p = 0.23 



1 M(P) H~ M(T,CFC) 



|M(E) [ 



p = 0.0002 



M(T,MacC)| 



|M(L) I- 



/? = 0.0002 



M(TJMRSA)| 



p = 0.015 



M(T,STAA) 



-t r~ 

6 9 
F value 



12 



— 1 

15 



Fig. 5. Representing the difference between 'type A' (thick boxes) 
and 4 type B* (thin boxes) models of the maximum specific growth 
rate by F-values. The respective P- values show the probability of 
that high F- value occurring by chance when the compared models 
are equal. 



p = 0.86 



< 

CD 

b 



M(T,CFC) | 



p = 8.4E-06 



p = 0.02 

M(T,MRSA)| 



H M(T,MacCj 



p = 0.02 



M(T,STAA) 



' i 1 — 

10 15 
F value 



i — 

20 



— I 

25 



Fig. 6. Representing the difference between the total (thick box) 
and 'type B' (thin boxes) models of the maximum specific growth 
rate by /"-values. 

test was carried out on the estimated growth rates 
(Fig. 6) as well as on the lag times (Fig. 7). 

Fig. 6 demonstrates that there was no significant 
difference between the model of the group P in the 
total mixture, M(T,CFC), and the total model, 
M(T,TSA). The growth of the rest of the groups was 
significantly slower than that of the total mixture. 
Hence, the group P was the dominant in the total 
mixture, and its model was the closest to M(T,TSA) 
while the model of the group E, M(T,MacC), was 
the furthest away. 



p = 0.66 



< 

tfi 

H 
H 



M(T,CFC) 



p = 0.006 



M(T,MacC) 



p-0.04 



M(T,MRSA)| 



p = 031 



M(T,STAA)| 



- I T~ 

2 4 
F value 



Fig. 7. Representing the difference between the total (thick box) 
and 'type B* (thin boxes) models of the lag by F-values. 
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Fig. 8. The generic model, M(G), fined on the specific growth 
rates of (i) the group P grown alone ('type A' model); (ii) the 
group P grown in the total mixture and selected by CFC agar 
('type B' model); (iii) the total mixture plated on TSA (total 
model). Legend as in Fig. la and Fig. lb. 



The obtained F-values for the lag, comparing 
'type B' models with the total model, are shown in 
Fig. 7. As the figure shows, the relation of the lag 
models are similar to that of the growth models. 

3.5. Creating a generic spoilage model for the 
maximum specific growth rate 

The generic spoilage model of the mixed popula- 
tion, M(G), based on the dominant group of the total 
mixture, was created by fitting those three datasets of 
In u,, for which the F test showed no significant 
differences (Fig. 8): 

• maximum specific growth rates of the group P, 
when grown in isolation; 

• maximum specific growth rates of the total mix- 
ture plated on CFC; 

• maximum specific growth rates of the total mix- 
ture plated on TSA. 

The obtained coefficients and their standard errors 
for this generic model arc shown in Table 5. 



4. Discussion 

Simple quantification and easy visualisation of the 



differences between predictive models is not new in 
the literature. The coefficients of predictive models, 
describing the growth rate of some Aspergillus spp. 
as a function of a v1 was used by Baranyi et al. 
(1997) as a measure of the relatedness of microbial 
species. Dengremont and Membre (1995), applied 
the percentage similarity between the specific growth 
rates of different Staphylococcus aureus strains to 
represent the distances between them. In our case, F 
test values were employed for a similar purpose. 

The F values quantifying, the discrepancy between 
'type A' and 'type B' models (Fig. 5) can establish 
an order between the groups in terms of level of 
inhibition. This order from the least to the most 
affected group can be summarised as group P< 
group B< group L< group E. 

In the mixed population (T), the closest model to 
M(T,TSA) is M(T,CFC), while M(T,MacC) is the 
most distant. Using the F values of Fig. 6 and Fig. 7 
in a similar manner as above, a rank can be 
established between the bacterial groups according to 
their contribution to the growth of the mixed popula- 
tion: group P> group B>group L> group E. 

Therefore, the dominant organism or organisms, in 
the mixed population, could be one or more of the 
strains of the group P, which were able to grow on 
CFC (P. putida, P. fragi, P. lundensis and /or S. 
putrefaciens). The rest of the spoilage groups were 
significantly inhibited in the mixed population. 

Dainty and Mackey (1992) raised the question 
whether the rapid growth of Pseudomonas spp. 
might be a sufficient reason for their dominance in 
some chilled food products. In our work, no signifi- 
cant differences were detected in the growth rates of 
the groups P, E and B when grown in isolation. Since 
the only difference in the growth conditions of the 
total mixture was the presence of the other groups, 
microbial interactions must be the cause of the 
reduction of the growth rates. We showed this 
reduction for all the organisms except for the group 
P, which was able to grow at the same rate in the 
total mixture as in isolation. 

Gill and Newton (1977) suggested that pseu- 
domonads may have a greater affinity for oxygen 
than facultatively anaerobic bacteria, which would 
result in the reduction of the growth rates of the 
latter group. However, this effect can be noticed at 
high concentration of pseudomonads only. 
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The ability of Pseudomonas species to transform 
glucose rapidly to a form not readily utilised by the 
competing organisms has been considered as a 
competitive advantage (Whiting et al. t 1976a,b). 
Indeed, the sequestering of glucose as gluconate 
might be one of the reasons for the success of 
Pseudomonas spp. (Nychas et a)., 1988). However, 
as Dainty and Mackey (1992) wrote, the inhibitory 
effect of the glucose transformation would be appar- 
ent only when glucose becomes limiting which is not 
the case in our broth system. 

Nychas et al. (1988) also mentioned that, beside 
the fast growth rates, the tolerance to acidity might 
be another reason for the success of Pseudomonas 
spp. However, Enterobacteriaceae, Moraxella, 
Aeromonas, Acinetobacter or Flavobacterium spp. 
are as tolerant to acidity as Pseudomonas spp. within 
the range pH 5.5-6.4 (Gill, 1986). 

Gram (1993) suggested that one of the causes for 
the bacterial selection in the microflora of some food 
products is the influence of microbial interactions 
such as competition or antagonism. She found that 
most strains of Shewanella putrefaciens were strong- 
ly inhibited by Pseudomonas strains isolated from 
spoiled and fresh fish. The inhibitory effect was 
associated with siderophore-mediated competition 
for iron. However, that inhibition was not observed 
in liquid media, only in agar assays, using cultures of 
pscudomonads pregrown to maximum cell numbers. 

Our work cannot clarify the reasons for the 
dominance of Pseudomonas spp. We can conclude 
only that the reason for this dominance is not 
exclusively their fast growth. Neither can it be that 
some substrates become limiting for the competing 
organisms since they were dominant, from the 
beginning of the growth of the total mixture, in a 
medium rich in nutrients, at relatively low bacterial 
concentrations. 

It is worth discussing the propagation of errors 
through the modelling procedure. In Table 3 and 
Table 4, the growth parameters were obtained 
through a non-linear curve-fitting procedure which 
also estimated the standard error of those parameters. 
These errors are expressed in the table as percentages 
of the estimated parameters themselves (relative 
errors). As can be seen, the relative errors for the 
maximum specific growth rates are fairly constant 
and are between 5% and 15%, dominantly below 



10%. The constant relative error justifies the use of 
the logarithm as a link function (McGullagh and 
Nelder, 1983) when the maximum specific growth 
rates were modelled as a function of temperature and 
pH. 

In Table 5, the standard errors of those fitted 
quadratic response surface models can be seen. 
Those errors are around 0.1 which indicates that the 
specific rates are estimated by about exp(0.1)=l.l 
factor accuracy. This means that the relative standard 
error of that estimation is about 10%, a result which 
is consistent with the overall picture of the relative 
errors of the specific rates given in Table 3 and Table 
4. In other words, the applied (in fact, purely 
empirical) quadratic polynomial model was 'good 
enough' in the studied region of the environmental 
factors because it did not introduce significant model 
error. 

The reason why we considered only the maximum 
specific growth rate in this comparative study, and 
not the lag time, is the doubtful reproducibility of the 
latter parameter (Shida et al., 1977; Buchanan and 
KJawitter, 1991; Hudson, 1993). 

Bacterial interactions and competition is exten- 
sively studied in literature. An inherent problem of 
these studies is that the quantity and quality of the 
available data do not match with the number of 
variables and the complexity of the mathematics 
required by mechanistic models. Panikov (1996) 
used four differential equations to model the steady- 
state growth and transient dynamics of a mixed 
population of pseudomonads and enterobacteria. To 
simulate the competition of 2-5 species for a 
common limiting substrate and the growth on several 
substrates, a model containing up to 22 differential 
equations was needed. The more accurately we want 
to describe a real system, the higher the number of 
necessary variables becomes, and the more complex 
the mathematical model will be. 

In this paper, our goal was not a mechanistic 
modelling of competition. Our approach was to 
quantify the concepts of 'dominance' and influence' 
of spoilage bacteria on each other. This aim was 
compatible with the data that we could measure with 
traditional microbiological techniques. However, fur- 
ther studies are necessary not only to identify the 
kind of interactions that can be established in a 
mixed culture, but also, to find a less laborious and 
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more accurate way of studying the behaviour of 
different organisms in heterogeneous populations. 
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